On the nature of ferromagnetism in Gai^Mn^As semiconductors 
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We have performed ab initio calculations within the density-functional theory for Gai-^Mn^As 
diluted semiconductors. Total energy results unambiguously show that a quasi-localized J, hole, with 
strong p-like character, surrounds the fully polarized Mn f d 5 -electrons. Calculations indicate that 
the holes form an essentially dispersionless impurity band, thus rendering effective-mass descriptions 
of hole states open to challenge. We obtain estimates both for the s = 1/2 hole and S = 5/2 Mn 
exchange coupling, and for the distance dependence of the effective Mn-Mn exchange interaction. 
Results demonstrate that the effective Mn-Mn coupling is always ferromagnetic, thus non-RKKY, 
and intermediated by the antiferromagnetic coupling of each Mn spin to the holes. 

PACS numbers: 71.55.Eq, 75.30.Hx, 75.50.Pp 



-a 
c 

o 
o 



> 

o 
m 
O 

-i— > 
c3 



i 

C 

O 

o 



X 



In the last few years, a considerable amount of work 
has been devoted to the study of diluted magnetic semi- 
conductors (DMS), since the possibility of manipulating 
both the charge and spin degrees of freedom of carri- 
ers in magnetic materials will change qualitatively the 
efficiency of spintronics devices. The discovery of hole- 
induced ferromagnetism in p-type (In,Mn)As systems 
was followed by the successful growth of ferromag- 
netic (Ga,Mn)As alloys Q. A quantitative understand- 
ing of the physics in these materials is therefore cru- 
cial, since ferromagnetic III-V alloys may be readily com- 
bined into semiconductor heterostruture systems, open- 
ing up a range of applications of optoelectronic devices 
through the combination of quantum and magnetic phe- 
nomena in these materials. However, in order to de- 
velop full-scale applications, one needs to elucidate sev- 
eral issues in relation to these systems. In the case of 
Gai-^Mn^As semiconductors, it is well known that Mn 
forms acceptors when in substitutional Ga lattice sites 
(MnGa)- One of the standing issues is related to the 
fact that the critical temperature and hole concentra- 
tion, as a function of Mn composition in Gai_ x Mn a As, 
are crucially dependent on the details of growth condi- 
tions [lIlMUHHSHmHHEi, even in as- 
grown samples; note, in particular, that for x ~ 5%, T c 
ranges from ~ 30 K to 110 K, with the highest tempera- 
ture Q of 110 K only being reproduced by other groups 
after post-growth annealing 0, 0, The measured 

hole concentration p as a function of x similarly reveals 
discrepancies both in trends and in order of magnitude; 
note that different experimental techniques (Hall resis- 
tance and Raman scattering 8J) applied to the same 
x = 0.083 sample yield values of p differing by as much 



as a factor of 10. 

From the theoretical point of view, different models 
have been proposed to describe the electronic and mag- 
netic properties of (Ga,Mn)As. The generally accepted 
view is that a given Mn ion interacts with the holes via 
a local antiferromagnetic Kondo-like exchange coupling 
Jhd = N/3 ■ (a 3 /4) (a is the GaAs lattice parameter) be- 
tween their magnetic moments [il llrjlll^ : this interaction 
is thought to lead to a polarization of the hole subsys- 
tem, which would then give rise to an effective ferromag- 
netic coupling between the Mn moments. Early attempts 
treated the holes within an sp parabolic-band effective- 
mass approximation; this approach is in conflict with re- 
cent photoemission studies 9] and infrared measurements 
[lij . which indicate that the holes form a relatively flat 
impurity band at the Fermi energy, instead of residing 
in an unaltered GaAs valence band. In particular, in- 
frared spectroscopic measurements estimated the hole ef- 
fective mass within the range 0.7 m e < m* < 15 m e for an 
x = 0.52 sample, and even larger values at higher dopings 
(lif . As a consequence, any treatment of holes within 
an effective-mass approximation, including attempts in 
the direction of incorporating aspects such as a Kohn- 
Luttinger treatment of the valence states [l^. l2(ij |. are 
quite clearly open to question. The same applies to es- 
timates of N/3 r~j 1.5 — 3 eV 0,^3; which are obtained 
from resistivity fits to experimental data resorting to a 
hole effective mass m* = 0.5 m e . 

Therefore, it is fundamental to have a highly accurate 
microscopic description of the electronic and magnetic 
properties of Gai-ajMn^As DMS. Here we perform a de- 
tailed ab initio study of the physical origin of the ob- 
served Mn-Mn ferromagnetic coupling. Our total energy 
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results provide unambiguous evidence that (i) the im- 
purity states are essentially localized and hence almost 
dispersionless, thus rendering any effective-mass picture 
inapplicable, and (ii) the effective coupling is always fer- 
romagnetic (thus non-RKKY), intermediated by an an- 
tiferromagnetic coupling of each Mn spin to the holes. 
In addition, we have also provided reliable values for the 
Mn-hole exchange coupling, as well as for the anisotropy- 
and direction-dependence of the effective (i.e., if the hole 
degrees of freedom were integrated out) coupling between 
Mn spins. 

We have performed total energy calculations based 
on the density-functional theory (DFT) within the 
generalized-gradient approximation (GGA) for the 
exchange-correlation potential, with the electron-ion in- 
teractions described using ultrasoft pseudopotentials [2lJ . 
A plane wave expansion up to 230 eV as implemented in 
the VASP code [23I was used, together with a 128-atom 
and a 250- atom fee supercell and the L-points for the 
Brillouin zone sampling. The positions of all atoms in 
the supercell were relaxed until all the forces components 
were smaller than 0.02 eV/A. We have also checked for 
spin-orbit effects through the projector augmented-wave 
(PAW) method |2^| and found that they may be safely 
neglected. 

Let us first consider the case of a single isolated MnQ a 
acceptor. From total energy calculations for neutral 
(MnQa) and negatively charged (MnQ a )~ defects, with a 
250-atom supercell, we obtain a localized acceptor level 
lying at 0.1 eV above the top of the valence band, in 
good agreement with the 0.11 eV experimental value |24j . 
Moreover, the ground state of the Mnc a defect is consis- 
tent with the picture of a J, hole interacting antiferromag- 
netically with the f S = 5/2 spin of the ^-configuration 
at the Mn site. The robustness of this state is illustrated 
by the fact that the ferromagnetic configuration with an 
| hole lies ~ 0.25 eV above the antiferromagnetic one. 
Figure^a) shows the difference m(r) = pt(r) — p±(r) for 
the MnQa defect, where p a is the total charge density in 
the cr-polarized channel. We note that near the MnQ a 
acceptor, the local magnetization has a strong f <i-like 
character, due essentially to the valence-band resonant 
d 5 electrons, whereas as one approaches its As neigh- 
bors, the character changes to j p-like. The signature of 
the above mentioned antiferromagnetic interaction con- 
sists of a sign change in m(r) as one moves from the Mn 
site to any of its neighboring As. In order to probe the 
extent of the electronic charge associated with the Mnc a 
defect, the difference between p-[ (r) + pi (r) calculated for 
the Mnca state and the GaAs host is depicted in Figure 
njb). One clearly notes the localized nature of the defect, 
as the charge difference is essentially confined within the 
region surrounding the Mn site. 

The origin of ferromagnetism in diluted Gai-^Mn^As 
semiconductors may be elucidated by focusing on inter- 
acting MnQa substitutional defects in a 128-atom super- 



cell, considering both a ferromagnetic as well as an anti- 
ferromagnetic allignment between the Mn spins. We have 
performed calculations for two Mn substitutional atoms 
in configurations corresponding to all inequivalent posi- 
tions within the supercell, i.e., Mn-Mn distances varying 
from 4.06 A up to 11.48 A. Our total energy results yield 
an unanbiguous Mn-Mn ferromagnetic ground state in all 
cases |25J. Figure El shows the net magnetization m(r) 
isosurfaces for two Mn defects in nearest-neighbor and 
next-nearest neighbor positions (for the Ga sublattice), 
with parallel and anti-parallel Mn spins. In Fig. |3 the 
Mn atoms are at the center of the spherical- like regions of 
m(r), whereas the p-like regions are always centered on 
As atoms. Note that, irrespective of the relative orienta- 
tion of the Mn-Mn spins, the antiferromagnetic coupling 
between the Mn and hole spins is always maintained. In 
the Mn-Mn antiferromagnetic case, this leads to the ap- 
pearence of nodes in m(r), which contributes to increase 
its energy relative to the ferromagnetic state. Another 
important feature illustrated in Fig. 0is the fact that, 
again, the Mn-Mn defect is essentially localized, although 
the magnetization density clearly spreads out from one 
Mnca site to the other. The picture that emerges is that 
of a cloud of I holes surrounding the substitutional Mn, 
with the distribution of quasi-localized holes giving rise 
to a nearly dispersionless impurity band, so that a de- 
scription via the effective-mass approximation would be 
inappropriate. 

From the difference between total energies of the ferro- 
and antiferromagnetic configurations of the s — 1/2 hole 
and the 5* = 5/2 Mn spins, assuming an interaction of the 
type N(3s ■ S, and taking into account the contribution 
due to images in neighboring supercells, one may cal- 
culate the strength of the Kondo-like antiferromagnetic 
exchange coupling N(3 ~ + 0.1 eV. This result contrasts 
with the previous estimates of N/3 ~ 1.5 — 3 eV |3. Il7|. 
We note that these larger values are open to question, as 
they are obtained via fittings, to experimental data, of 
results based on the effective-mass approximation with a 
hole effective mass m* = 0.5 m e . We have also estimated, 
from total energy calculations, the effective exchange 
coupling between pairs of S = 5/2 Mn spins, J Mn-Mn, 
as a function of the Mn-Mn distance, for all inequivalent 
pair positions within the supercell (for a Mn-Mn interac- 
tion modelled through J Mn-Mn Sm^ • S M n 3 )- In Fig. 
we display the theoretical predictions for J Mn-Mn- The 
results clearly show that the coupling between the Mn 
spins is always ferromagnetic, irrespective of their rel- 
ative distance. As it is well known that the bare cou- 
pling between two Mn spins should be antiferromagnetic, 
one concludes that the resulting Mn-Mn ferromagnetic 
effective coupling, in Gai-^Mn^As, is essentially inter- 
mediated by the antiferromagnetic coupling of each Mn 
spin to the quasi-localized holes. Also, the observed non- 
monotonic behavior of JMn—Mn{ r ) should be attributed 
to the anisotropic character of the effective interaction, 



3 



as it may be inferred from Fig. |2J The inescapable con- 
clusion from Fig. [3|is that an RKKY description for the 
Mn-Mn interaction in Gai_ x Mn.jAs is ruled out. Here we 
note that this result contrasts with the RKKY coupling 
obtained by Zhao et al. 26] in the case of Mn^Gei-j. 

The data in Fig.|3|can be used to estimate the critical 
temperature at a given Mn concentration, which, within 
a simple mean-field theory is given by 
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where Ji, J2,... stand for first-, second-,..., -neighbor 
interactions, the values of which are given in Fig. [31 and 
the ~zi are the corresponding configurationally-averaged 
coordination numbers in the Ga sublattice. If we take 
z7 = fzi, with / being the fraction of active Mn sites, 
which is roughly x/3, one obtains T c — 4 x 10 3 a; K, in 
qualitative agreement with the low-density behavior ob- 
served experimentally 0]. 

One should mention that the discrepancies between 
different experimental data indicate that the mag- 
netic, structural, and electronic properties of as-grown 
Gai-a-Mn^As alloys are extremely sensitive to the ac- 
tual molecular-beam epitaxy (MBE) growth conditions, 
such as, for example, growth temperature and beam flux 
ratios. In fact, this is to be expected, since the holes 
provided by the Mnc a acceptors may be compensated 
by defects such as arsenic antisites (Asc a ) donors, Mn 
interstitials (Mn/), Mnj-MnQa pairs, MnAs complexes, 
etc. Of course, the situation with respect to annealed 
Gai-xMn^As samples is rather more complicated, and a 
proper analysis of experimental measurements performed 
in after-annealing samples must involve a realistic mod- 
elling of diffusion processes involving several defects, for- 
mation of random precipitates, clustering effects, etc. 
Therefore, we emphasize that a proper understanding of 
the physics of Gai-^Mn^As alloys must involve a micro- 
scopic description of the effect of different defects on their 
electronic and magnetic properties. 

In conclusion, we have provided a detailed ab initio 
study of the physical origin of the Mn-Mn ferromagnetic 
coupling, by considering isolated MnQ a defects, as well 
as two substitutional Mn per supercell, in various rela- 
tive positions. Our total energy calculations provide un- 
ambiguous evidence that the effective coupling is always 
ferromagnetic, and thus non-RKKY, and intermediated 
by an antiferromagnetic coupling of each Mn spin to the 
quasi-localized holes. 
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Figure 1: Isosurfaces for (a) the net local magnetization 
m(r) = P|(r) — pj.(r) for the MnGa defect, and (b) the differ- 
ence between P|(r) + pi(r) calculated for the MnGa ground 
state and the GaAs host. The green surface corresponds to a 
value of + 0.004 e/A 3 , and the blue surface to - 0.004 e/A 3 ; 
for comparison, note that a uniform charge density of 1 elec- 
tron/unit cell in GaAs corresponds to + 0.022 e/A , with e 
being the electron charge. The black (red) spheres denote the 
Ga (As) atoms. 



Figure 3: Exchange coupling Jmu-Mu between Mn atoms in 
Gai-ajMn^As alloys, versus the Mn-Mn distance. Open dots 
are our calculated results whereas the full line is a guide to 
the eye. 



Figure 2: Isosurfaces for the net local magnetization m(r) = 
Pl (r) — p| (r) in the case of two MnGa defects; the color code 
and isosurface values are the same as in Fig. 1. In (a) and (b) 
[(c) and (d)] the two Mn are nearest neighbors [next-nearest 
neighbors] with their S = 5/2 spins alligned parallel and anti- 
parallel, respectively. 
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